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The distribution of entangled quantum states of light 
over long distances is a major challenge in the field of 
quantum information. Optical losses, phase diffusion and 
mixing with thermal states lead to decoherence and de- 
stroy the nonclassical states after some finite transmis- 
sion line length. Quantum repeater protocolsi^ combin- 
ing quantum memory^, entanglement distillation^ and 
entanglement swapping* were proposed to overcome this 
problem. Here we report on the first experimental demon- 
stration of entanglement distillation in the continuous 
variable regime' -'. Entangled squeezed states were first 
disturbed by random phase fluctuations and then distflled 
and purified using interference on beam splitters and ho- 
modyne detection. Measurements of covariance matri- 
ces clearly indicated a regained strength of entanglement 
and purity of the distiUed state. Contrasting previous 
demonstrations in the complementary discrete variable 
regimeiSili, our scheme achieved the actual preparation 
of the distUled states, which might therefore be used to 
improve the quaUty of down-stream applications such as 
quantum teleportation^ 

Quantum information makes use of the special properties 
of quantum states in order to improve the quality of com- 
munication and information processing tasks. Generally, a 
quantized field can be described by the number operator or 
alternatively by two non-commuting position and momentum- 
like operators. The corresponding measurement results have 
either discrete or continuous spectra and form the basis of 
discrete variable (DV) or continuous variable (CV) quantum 
information, respectively. In the regime of continuous vari- 
ables, entangled states of light can be deterministically gen- 
erated in optical parametric amplifiers (OPAs), precisely ma- 
nipulated with linear optics, and measured with very high effi- 
ciency in balanced homodyne detectors. These entangled two- 
mode squeezed states show Gaussian probability distributions 
and were utilized for quantum teleportatior^lS and entangle- 
ment swappingi^iii. Entangled states of the collective spins of 
two atomic ensembles analogous to two-mode squeezed states 
have been generated^^, storage of quantum states of light in 
an atomic memory has been demonstrated^ and teleporta- 
tion from light onto an atomic ensemble has been reported'^. 
High-speed quantum cryptography with coherent light beams 
and homodyne detection has been demonstrated'^. All these 
spectacular achievements reveal the great potential of this ap- 
proach to quantum information processing. 

A missing piece in this toolbox has been a feasible pro- 
tocol for entanglement distillation and purification. Entan- 
glement distillation"^'^ extracts from several shared copies of 
weakly entangled mixed states a single copy of a highly en- 
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FIG. 1: Experimental setup. OPA: optical parametric amplifier, 
BSd: distillation beam splitter, BSe: entangling beam splitter, BHD: 
balanced homodyne detector. (Jli indicate independently fluctuating 
optical path lengths. The combination of both BHDt.a/b provided 
the trigger signal for successful distillation, whereas BHDv,a/b 
were used for its independent verification. The latter were not neces- 
sary for the distillation protocol to work. 



tangled state using only local quantum operations and classi- 
cal communication between the two parties sharing the states. 
This turned out to be a very challenging task for CV states 
since it was proved that it is impossible to distill Gaussian 
entangled states by means of the experimentally accessible 
Gaussian operation s'^-'^ . However, a whole class of impor- 
tant decoherence processes give rise to non-Gaussian noise 
and therefore produce non-Gaussian entangled states. It has 
been shown^ that in this case the entanglement distillation can 
be carried out using only interference on beam splitters, bal- 
anced homodyne detection and conditioning on the measure- 
ment outcomes. This was recently experimentally confirmed 
by successful demonstrations of distillation and purification 
protocols for single squeezed modes that suffered from de- 
Gaussifying noise 
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This work experimentally demonstrates the first distilla- 
tion protocol for CV entangled states. Our protocol uses 
two copies of phase-diffused22i2I, therefore mixed, entangled 
states exhibiting a positive, but non-Gaussian Wigner func- 
tion. Our protocol enhances the entanglement and purity of 
the decohered states and represents a single step of an itera- 
tive Gaussification schem^^ that asymptotically converts any 
input state into a Gaussian one. Moreover, if combined with a 
single de-Gaussifying operation such as the recently demon- 
strated single-photon subtraction from squeezed beamssi^^i^^, 
it would provide a generic continuous variable entanglement 
purification and distillation scheme--^, that is capable, for in- 
stance, of suppressing the detrimental effect of losses in quan- 
tum state transmission. 

A schematic sketch of our protocol is shown in Figure [T| 
Two optical parametric amplifiers (OPAs) provided two con- 
tinuous wave light fields that carried squeezed states of light. 
Both squeezed states were mixed with vacuum states on beam 
splitters with 50% power reflectivity in order to prepare two 
copies of so-called v-class2i entangled states. Entanglement 
prepared in this way is not the strongest possible but is effec- 
tual for this proof of principle experiment. All four resulting 
beams of the two copies of entangled states were transmit- 
ted to two parties Alice (A) and Bob (B) through four chan- 
nels exhibiting independent phase noise. The noisy channels 
were realized by quasi-random electro-mechanical actuation 
of mirror positions in the beam paths22i21 in order to mimic 
the phase noise introduced for example in optical fibres. The 
phase fluctuations applied exhibited a Gaussian distribution, 
hence the standard deviation of the noise provided a com- 
plete characterization of its strength. 

Alice and Bob each received two beams which they over- 
lapped on a balanced beam splitter (BSd.a and BSd,b in 
Fig. [Til. The mean phase was controlled such that the initial 
quadrature phases lined up. The two output ports of each 
beam splitter were connected to a total of four balanced ho- 
modyne detectors (BHDs). Each of the BHDs could be set 
to observe an arbitrary quadrature, in particular the amplitude 
quadrature X or the conjugate phase quadrature P. The de- 
tectors BHDt.a and BHDt.b were used to generate trigger 
events discriminating the success or failure of preparation a 
disitlled state, while BHDv,a and BHDv,b were used to in- 
dependently verify the performance of the protocol. The latter 
two were therefore not part of the actual protocol and could 
be replaced by an arbitrary experiment. The fringe visibili- 
ties on all four BHDs as well as on (BSd,a/b) were between 
97.1% and 98.2%. In our experiment both trigger detectors 
BHDt A/B were set to measure the initially squeezed ampli- 
tude quadratures Xt,a and Xt.b, respectively. We chose the 
trigger condition determining the success of the distillation 
protocol in the form^ 



l^T.A + Xtr\ < 



(1) 



where Q was a certain threshold whose value could be tuned 
in order to vary the selectivity of the protocol. 

Figure |2^) shows the variance of the nonlocal EPR-like 
quadrature operator {Xv,a + Xv.b) of the successfully dis- 
tilled state and its corresponding probability of success- 
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FIG. 2: Nonlocal quadrature variance and total variance of the 
distilled states, (a) Variance of (Xv.a + Xv,b) and corresponding 
success rates versus threshold value Q for different strengths of the 
phase fluctuations apn- o : measurement, - : simulation, (b) The total 
variance T plotted against the success rate can be seen as measure for 
the performance of our protocol. The grey lines indicate the vacuum 
noise level or separability boundary, respectively. For a success rate 
as high as 0.5 the protocol deploys nearly its full potential. 



ful preparation versus threshold value Q for three differ- 
ent strengths of the phase fluctuations apn- The lower the 
threshold was set, the more selectively our protocol worked. 
The variance of the orthogonal EPR-like quadrature operator 
(Pv.A — Pv.b) was also reduced (not shown). In this case 
the effect was rather small and limited by the vacuum noise 
value of 1/2, because v-clas922. entangled states were used in 
our experiment. 

Figure I2J)) shows the total variance T of the distilled state 
versus its preparation success rate for five different strengths 
of phase noise. The total variance quantifies the quantum cor- 
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FIG. 3: Reconstructed covariance matrices of distilled entangled 
states. Displayed is the lower left triangle of this matrix containing 
its ten significant elements. The sub-bars represent the results of our 
distillation protocol for different threshold values Q in the case of 
o"pn = 0.497. The base area of each sub-bar represents the corre- 
sponding success rate. The vacuum reference for this plot is given 
by the unit matrix. 



relations between two modes and the degree of inseparability 
of the joint state (see Methods). If X < 1 then the state is not 
separable and therefore entangled^^. We can see that the dis- 
tillation and purification protocol reduces the total variance 
T of the phase-diffused state. For the two strongest levels 
of phase noise, the total variance shows no nonclassical be- 
haviour without the distillation protocol, but was reduced be- 
low the unity boundary after applying the distillation protocol. 
Lower values of Q result in a stronger distillation effect but 
also in a reduced success probability. It is promising to find 
that for a success rate as high as 0.5 the protocol had nearly 
developed its full potential. Note that the total variance of the 
entangled state before applying phase noise was evaluated to 
beX = 0.725. 

Since the homodyne detectors can measure arbitrary 
quadrature phases, we performed a tomographic reconstruc- 
tion of the covariance matrices of the distilled two-mode en- 
tangled states. We followed the reconstruction procedure that 
is described in detail in our previous worb21 and is based on 
joint measurements of quadrature operators of modes V,A and 
V,B for four different settings of the BHDs. From this data 
we determined the most significant eight of the ten indepen- 
dent parameters of the covariance matrix, namely the vari- 
ances of four quadratures Xva, Pv,a, Xyb, Py,b and covari- 
ances between all pairs of quadratures of Alice and Bob. The 
intra-modal correlations were neglected, because the devia- 
tion from zero was too small to be measured within the given 
accuracy of the quadrature phase. The resulting covariance 
matrices are plotted in Fig. [3] for ten different success rates. 
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FIG. 4: Determinant of covariance matrix of distilled state The 

determinant D is plotted versus the success rate of the distillation 
protocols for five different levels of phase noise. The decrease of the 
determinant D towards the left is an indication of increased purity 
and Gaussification of the distilled state. 



With decreasing success rate the distillation becomes stronger 
as witnessed by the reduction of the quadrature variances, 
i.e. diagonal elements of the covariance matrix. Moreover, 
the anti-correlation between Xv a and Xvb was slightly en- 
hanced. Consequently, the squeezing of the nonlocal quadra- 
ture Xv,A + ^v,B was enhanced by the distillation. 

The determinant D = det(7) is an important characteristic 
of a covariance matrix 7. The purity V of a Gaussian state 
p defined as V = Tr[p2] is given by "P = 1//D. In Fig.g] 
we plot the dependence of the determinant D on the success 
rate of the distillation protocol. The distillation reduces D 
which is a very strong indication of the increased purity of the 
(still slightly non-Gaussian) distilled state, because a general- 
ized Heisenberg uncertainty relation implies that D > 1 and 
D = 1 holds only for pure Gaussian states. The anticipation 
that the distilled states should not only show purification but 
also Gaussification was indeed confirmed by fitting Gaussian 
functions to the measurement data before and after the distil- 
lation protocol. 

The protocol demonstrated here is able to distill entangled 
states from a decohered ensemble that suffered from phase 
noise in the transmission channel e.g. due to thermally excited 
refractive index fluctuations or Brillouin scattering. Because 
of its quantum nature, our protocol can counteract phase diffu- 
sion in cases where any classical channel probing would fail, 
for example if the phase fluctuations are intensity dependent 
and arise from quantum fluctuations in the photon number of 
the transmitted states. Our protocol provides two open ports 
that output the distilled entangled states, and is therefore un- 
conditionally useful for arbitrary down-stream quantum infor- 
mation applications that involve the second order moments of 
quadrature operators. An example is a teleportation protocol 
that teleports Gaussian statesi2,. The reported entanglement 
distillation, purification and Gaussification protocol can be 
iterated^'** and combined with already experimentally demon- 
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strated single-photon subtraction2i2^i2^, quantum memorji^ 
and entanglement swappin g'^-'"^ to build a continuous-variable 
quantum repeater. Our experiment is thus an important en- 
abling step towards truly long-distance quantum communica- 
tion with continuous variables. 



METHODS 

Generation and detection of squeezed states Both opti- 
cal parametric amplifiers were constructed from second order 
nonlinear crystals (MgO:LiNb03) inside a degenerate doubly 
resonant cavity^l and were pumped with frequency doubled 
laser beams from a monolithic solid state laser (Nd:YAG) 
operated at 1064 nm. The OPAs emitted continuous wave 
squeezed beams with about 4.5 dB of squeezing and 8 dB anti- 
squeezing at Fourier frequencies around 7 MHz. 

All four output beams of the distillation protocol were mea- 
sured by balanced homodyne detectors (BHDs) each consist- 
ing of a balanced beam splitter, a coherent local oscillator 
beam at 1064 nm and two photodiodes. The photocurrents 
produced by the photodiodes were subtracted and amplified 
by home-made electronics. The signals from the four BHDs 
were electronically mixed with a 7 MHz local oscillator The 
intermediate signals were then anti-alias filtered at 400 kHz, 
synchronously sampled with 1 MHz and post-processed to ac- 
complish and verify the distillation protocol. 

Entanglement detection and quantification The entan- 
glement of the distilled two-mode squeezed output state can 
be conveniently certified by the Duan criterion^^. This cri- 



terion involves the sum of variances of two nonlocal, EPR- 
like quadrature operators, the so-called total variance, X — 
((AX+,v)') + ((AP_,v)'), where X+ 
and -P_,v = ^v.A ^ A'.b- The state is entangled if X < 1, 
where the vacuum state corresponds to a value of X = 1. 
The total variance T provides a quantitative measure of the 
quantum correlations between the quadratures of the distilled 
states. For symmetric two-mode squeezed Gaussian states I 
fully quantifies the entanglement of the state-''. It has been 
shown numerically- that even for the non-Gaussian phase dif- 
fused squeezed states there is a one-to-one correspondence 
between the reduction of the total variance and the increase 
of the logarithmic negativity, which is a computable entangle- 
ment measure-'". 
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